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ABSTRACT 

High energy gamma-ray emission from two nearby bright starburst galaxies, M82 and NGC 253, have re- 
cently been detected by Fermi, H.E.S.S., and VERITAS. Since starburst galaxies have a high star formation 
rate and plenty of dust in the central starburst region, infrared emissions are strong there. Gamma-ray photons 
are absorbed by the interstellar radiation field photons via electron and positron pair creation. The generated 
electron and positron pairs up scatter the interstellar photons to very high energy gamma-ray photons via cas- 
cade emission through inverse Compton scattering. In this paper, we evaluate the contribution of this cascade 
emission to the gamma-ray spectra of M82 and NGC 253. Although it would be difficult to see direct gamma- 
ray evidence of cosmic-rays with an energy > 10 TeV due to the gamma-ray attenuation, the resulting cascade 
emission would be indirect evidence. By including the cascade component, we find that the total flux above 
1 TeV increases ~ 18% and ^45% compared with the absorbed flux assuming the maximum kinetic proton 
energy as 45.3 TeV and 512 TeV, respectively. Future gamma-ray observatories such as CTA would be able 
to see the indirect evidence of cosmic-ray with an energy > 10 TeV by comparing with theoretical emission 
models including this cascade effect. 

Subject headings: cosmic rays - galaxies: individual (M82, NGC 253) - gamma rays: general 



1. INTRODUCTION 

Supernova remn ant (SNR) shocks are believed to accel- 
erate cosmic-rays dGinzburg & Svrovatskiilll964T ; lHavakawal 
1969). Although this was confirmed by the detection of non- 
therm al emissions from the Galactic SNRs in X -ray and 7- 
rays dKovama et all 1 1995b lAharonian et alj|2004l) . the cause 
of the highest energy galactic cosmic -rays is still under dis- 
cussion. The highest energy of the Galactic cosmic -ray are 
thought to be comparable wi th the knee energy of the c osmic - 
ray spectrum, ~ 10 15 5 eV dNagano & Watsod 1200 0). Very 
high energy (VHE; >30 GeV) 7-ray observations of SNRs 
shows that SNRs have a cut-off around 10 TeV in the 7- 
ray spectra or ve ry soft injection proton spectra (see e.g. 
lAbdo et aT]|2010bh . Thus, it is not possible to simply expect 
that the SNR shocks accelerate cosmic -rays up to the energy 
around the knee energy region. Further investigation about 
these problems is currently required. 

Since starburst galaxies have a high star formation rate 
(i.e. high supernova rate) and plenty of gas, they are a 
candidate for proving the cosmic -ray acceleration scenario 
as well as the Galactic SNRs. Many papers have studied the 
7-ray emission mechanism in starburst galaxies dvbelk et al] 
19891: lAkvuz et alJl!99lb iPaelione et al] [ 19961: iTorresI ' 
Domingo-Santamana & Torres 2005: iPersic et all 
de Cea del Pozo et alJ 120091: iRephaeli et all 120101) Very 
recently, high energy 7-ray emissions from the nearby bright 
starburst galaxies, M82 and NGC 253, have been detected 
by Fermi H.E.S S. and VERITAS dAbdoetalj 1201 Oat 
lAcero et all 120091: IVERITAS Collaboration etaiT l2009l) . 
These detections support the idea of 7-rays generated by 
cosmic-ray propagation in starburst galaxies. Observation- 
ally, however, the existence of high energy emissions > 10 
TeV has not been satisfactorily constrained yet, although it 
would be a probe of the existence of cosmic -rays with an 
energy > 10 TeV. 

It is well known that VHE 7-rays propagating through the 
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universe are absorbed by the cosmic optical-infrared back- 
grou nd radiation via positron and electron, e + e~, pair produc- 
tion ( Gould &Schre der 1966[ |jeiievl[T966h . This 7-ray atten- 
uation process also occurs i n starburst galaxies which have 
plenty of infrared photons dDomingo-Santamarfa & Torres! 
I2005tlde Cea del Pozo et alJ2009l) . It might be difficult to de- 
tect > 10 TeV emissions from starburst galaxies due to this 
attenuation process. 

Here, created electron-positron pairs by attenuation up- 
scatter the interstellar photons to VHE 7-ray photons 
via the inverse Compton scatteri ng process. This is 
the so cal l ed ca s cade emi s sion (lAharonian et alJ | 1994t 
IWang et al] 120011: iDai et alJ 120021; iRazzaaue eTall 2004 
| Anddl2004lMurase et al.] 2007; Kneis ke & MannheirrJl2008l: 
llnoue & Totanill2009l:IVentersll2010h . If the energy of the in- 
trinsic 7-rays extends to >10 TeV, such a cascade emission 
would appear in the VHE band. This cascade emission would 
be indirect evidence of the existence of intrinsic > 10 TeV 
emissions. Therefore, we need to take into account cascade 
emissions to investigate the existence of cosmic -ray with an 
energy > 10 TeV. In the past studies on starburst galaxies, 
however, the effect of cascade emission was not taken in to 
account. 

In this paper, we estimate the contribution of the cascade 
emission to the 7-ray spectra for M82 and NGC 523. To in- 
vestigate this effect, we calculate the electron-positron pair 
creation optical depths in these starburst galaxies. We also 
compare the expected spectra including absorption and cas- 
cade effects with observed 7-ray spectra. 

2. ABSORPTION AND CASCADE EMISSION MODELING 

It is well known that VHE photons from high redshifts 
are absorbed by the interaction with the intergalactic optical- 
infrared radiation field via electron -positron pair creation 
dGould & Schrederll 19661: 1 Jellevlfl 9661) . It is also theoretically 
expected that these created pairs would scatter the cosmic 
microwave background radiation as secondary 7-ray emis- 
sions (the so called cascade emission) in the case of blazars, 
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20021: Razzaqu e et al.l | 2004t lAnddl2004t iMurase et al.ll200 
Kneiske & Mannheim! l2008h llnoue & Totanil 120091: IVenter; 
20101) . We investigate this cascade emission effect in two 
nearby starburst galaxies, M82 and NGC 253. 

2.1. Absorption ofVHE j-ray photons in starburst galaxies 

Very recently two nearby bright starburst galaxies, M82 
and NGC 253, have been detected with 7-ray observations 
by Fermi H.E.S. S. and VERITAS dAbdo et all 1201 Oat 
lAcero et all 120091: I VERITAS Collaboration et all 120091 
These 7-ray detected starburst galaxies would not suffer from 
7-ray absorption with the intergalactic radiation field, since 
they are very clos e to us. The distance to M82 and NGC 253 
is 3.6 ±0.3 Mpc dFreedmanetal.il 19941) and 3.9 ±0.4 Mpc 
dKarachentsev et al.ll2003l) . respectively. The starburst region 
is located in the inner part of the galaxy with a s ize of 500 pc 
and 280 pc for M82 a nd NGC 253, respectively dMavva et all 
l2006t IUlvestadll2000l) . 

The central starburst region in these galaxies where SNRs 
are produced is theoretically expected to emit most of the 7- 
ray photons (Domingo-Santamaria & Torres 2005, hereinafter 
DST05; de Cea del Pozo et al. 2009, hereinafter dCdP09). In 
the case of the Large Magellanic Cloud, the association be- 
tween starburst region and 7- ray emission region is observa- 
tionally confirmed by Fermi dAbdo et al.l [20 10c). It is also 
known that there are plenty of interstellar radiation field pho- 
tons from star forming activities. Therefore, we need to take 
into account the 7-ray absorption effect in the central star for- 
m ation region. 

ISiebenmorgen & KriigeH (120071 hereinafter SK07) devel- 
oped optical-infrared spectral energy distribution (SED) mod- 
els for the nuclei of starburst galaxies and ultra luminous 
infrared galaxies by computing a radiative transfer SED of 
spherical, dusty galactic nuclei. They include silicates, amor- 
phous carbon, graphite grains, and polycyclic aromatic hydro- 
carbons in the form of dust. Fig. Q]shows the optical-infrared 
spectral energy distribution (SED) of M82 and NGC 253. We 
adopt their spectrum model for M82 and NGC 253 shown in 
Fig. 2 and 3 in their paper, which is consistent with the data. 
In the case of M82, we add a blackb ody component (T = 2500 
K) to fit the data below 5 /xm as in ISiebenmorgen & Krugell 
d2007l) . which was not taken into account in previous works. 
For comparison, we also show the interstellar SED model 
in dCdP09 and DST05 for M82 and NGC 253, respectively. 
Since SK07 models well reproduce the observed data over 
a wide range of wavelengths for both galaxies, hereafter we 
use the SK07 optical-infrared interstellar SED models as our 
baseline nuclei SED models. 

Fig. |2] shows the expected optical depth, t 77 , of M82 
for SK07 and dCdP09 models and NGC 253 for SK07 and 
DST05 models. We do not take into account the interaction of 
a 7-ray wi th a nuclei, since this would not absor b 7-rays sig- 
nificantly (Domingo-Santamaria & Torres 2005). The 7-7 
absorption becomes significant (i.e. r 77 = 1) above ~ 5 TeV 



and ~ 9 TeV for M82 and NGC 253, respectively. 
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Persic et al.ll2008l:lde Cea del Pozo et alJl2009l:lRephaeli et alj 



I20T0I) . In this paper, we substitute the 7ro-meson decay 
spectrum for the intrin sic 7-ray spectrum, dN /dE^, us- 
ing code provided by Karlsson & Kamael d2008l) (see also 
iKamae et al.ll2005l 120 06). For the total inclusive inelastic p-p 
cross section, they include the non-diffractive (with Feymann 
scaling violation) and diffractive components, plus the 
A(1232) and Res(1600) resonance excitation contributions. 
The normalization of the total 7-ray flux including absorption 
and cascade effects is adjusted to the observed Fermi photon 
flux data in the 0.1-5 GeV band. The intrinsic 7-ray spectrum 
is given by the highest proton kinetic energy, T f 



p.max? 



in 



units of [TeV] and the spectral index of the intrinsic proton 
spectrum, T, where we set dN p /d~f p oc 7~ r , dN p /dj p the 
proton flux, 7^ the Lorentz factor of protons. We do not take 
into account inverse Compton scattering and bremsstrahlung 
radiation, since their contribution is expected to be less than 
10 % above 0.1 GeV (DST05, dCdP09). Although we need 
to solve the cosmic-ray propagation in starburst galaxies for 
detailed calculations, we adopt a simple analytical model to 
see the absorption and cascade effects. 

For the calculation of the cascade emission, we ignore mag- 
netic fields, since the 7-ray emissions from starburst galaxies 
are isotropic and continuous. We calculate the resulting cas- 
cade emissions as follows: 



dN c 
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Ando]|20o4lMurase et al.ll20q7t iKneiske & MannheirrJ l2008: 
Inoue & Totanill2009t IVentersll2010h . where f/c is the inverse 
Compton (IC) cooling time scale and the electron injection 
spectrum is 
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and the scattered photon spectrum per unit time by IC scat- 
tering is: 
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with/(;t) = 2x\n{x)+x+\-2x 2 , (0 <x < 1) andx = E- f . i /4^e. 
Here, £ 7 , = 2j e m e c 2 is the energy of intrinsic photons, 
dNj/dEyj is the intrinsic 7-ray spectrum, ro the classical 
electron radius, dn/de(e) is the photon density at the star- 
burst region. The integration region over the Lorentz fac- 
tor, 7 e , is 7 e . min < 7 e < 7 e , max , 7<-.max = £' 7 ,m a x/2 and 7 f , !min = 

max[m e c 2 /2e,(£ 7 /e) 1 ' /2 /2], where £ 7l max is the maximum en- 
ergy of 7-ray photons. 

We iteratively calculate Eq. [T|by substituting dN c /dE^(l - 
g-rvr) m or( j er to include IC scattering due to generated pairs 
from reabsorbed secondary photons. 

3. RESULTS 

Fig. [3] shows the 7-ray spectra of M82 and NGC 253 
in units of [MeV/cm 2 /s]. We show the intrinsic (the spec- 
trum without taking into account absorption), absorbed, and 
cascade components as well as total (absorbed + cascade) 
spectra. For the intrinsic spectra, we adopt (r,7J, max ) = 
(2.4,512 TeV) and (2.3,512 TeV) for M82 and NGC 253 
as our standard models, respectively. In the case of M82, 
our 7-ray spectrum nicely fits to the observed data in the 
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FIG. 1 .— Left panel: SED of M82. The solid and dashed c urve represents [Siebenmorgen & Kriigel ( 2007 a SK07 ) andlde C ea del Pozo et al] 120091 
dCdP09) SED model, respectively. Data references (near infrared: Johnson ( 1966); Kleinmann & Low ( 1970); Aaronson ( 1977); Jarrett et al. 12003); far in- 
frared: IRieke & Lowlfi97 2); Telesco & Harper 1 1980); Rieke et al. 1 1980); Klein et al. ( 1988); Telesco & Gezari ( 1992); Forster Schr eiber et aTM 2003); submm: 
I Jaffeetalj U9 84); Krugel et al. (1990); Hughes eTall <1990L I1994IV ). Right panel: SED of NGC 253 The solid and dashed curve represents SK 07 and 
Domingo-Santamaria & Torres 12005, DST05) SED model, respectively. Data references (near infrared: Rieke & Low 1 1975); Radovich et al. (2001); far in- 
frared: [Telesco & Harper ( 1980); Melo et al. 1 2002]); submm: IRieke et al] 1 1973): Hildebran d et al] fi977 ): Elias et al. 11978); Chini et al. H1984ft ). 
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FIG. 2. — The optical depth of pair production, r 77 , in starburst galaxies 
as a function of the photon energy. Solid, dashed, dotted and dot-dashed 
curve corresponds to t 77 in SK07 model for M82, dCdP09 model for M82, 
SK07 model for NGC 253 and DST05 model for NGC 253, respectively. The 
horizontal line marks the level of optical depth t 11 = 1 . 



range of the uncertainty of the data by taking into account 
the 7-7 absorption and the cascade emissions. However, 
in the case of NGC 253, it is difficult to explain the Fermi 
data point at 10 GeV, even if we take into account the sec- 
ondary 7-ray populations. Since the absorption is effective 
above 10 TeV and the seed photon density has a peak around 
0.01 eV, the typical energy of the secondary 7-ray, E c , is 
E c ~ (10TeV/2m e c 2 ) 2 x O.OleV - ITeV. Thus, cascade emis- 
sion do not contribute at ~10 GeV, rather, closer to around 
~ 1 TeV. By including cascade emissions, the total 7-ray flux 
above 1 TeV becomes 43% and 46% greater than the absorbed 
flux for M82 and NGC 253, respectively. 

Fig. @]also shows the expected 7-ray spectra of these two 
starburst galaxies but changing the intrinsic spectral index of 
proton, r, the highest proton kinetic energy, T p , max , and in- 
terstellar radiation field SED models for comparison. In the 
case of M82, the total 7-ray flux above 1 TeV becomes 48 
% and 22 % greater than the absorbed flux with (T, T p max ) = 
(2.3, 512 TeV) and (2.4,45.3 TeV) model, respectively. When 
we use the dCdP09 interstellar radiation field SED model with 
the standard parameters, the total flux above 1 TeV is 47 % 
greater than the absorbed flux. In the case of NGC 253, the 



total 7-ray flux above 1 TeV becomes 52 % and 14 % greater 
than the absorbed flux with (r,r pmax ) = (2.2,512 TeV) and 
(2.3,45.3 TeV) model, respectively. When we use the DST05 
SED model with the standard parameters, the total flux above 
1 TeV is 37 % greater than the absorbed flux. Therefore, even 
if we use other interstellar radiation field model which do not 
fit to the near infrared observed data, the effect of the cascade 
emission would not be changed significantly. This is because 
the near infrared photon density is not high enough to absorb 
VHE 7-rays effectively. 

It would be difficult to explain the data by the standard pa- 
rameters, but (r, T p , max ) = (2.2,4.0 TeV) model would be able 
to explain the data. We should note, however, note that, in this 
case, cascade emissions for this model would not contribute to 
the VHE spectrum since >TeV emission is weak. However, 
Tp.max = 4.0 TeV model for M82 is inconsistent with the data. 
We will need to wait for much longer time integration data 
from Fermi and future TeV 7-ray observatories to discuss this 
point in greater detail. 

4. DISCUSSION 
4.1. Implications for future VHE observations 

To see the feature of the cascade emission, it will be im- 
portant to detect these two nearby starburst galaxies with a 
high signal-to-noise ratio with future VHE observations. The 
Cherenkov Telescope Array (CTA) is planned for the next 
generation imaging atmospheric Cherenkov telescope (IACT) 
1 . The sensitivity and energy range of CTA will be improved 
by one order of magnitude compared to current IACTs such 
as H.E.S.S, MAGIC, and VERITAS. Here we investigate the 
required observational time to see the absorption and cascade 
signature with high significance by fut ure CTA observati ons 
following the argument in §3 and §4 in Ino ue et all d2010l) . 

First, the required observational time to detect the integral 
7-ray flux above 1 TeV with 5 a is 1.3 hrs for M82 and 1.0 
hr for NGC 253, respectively. Next, to see the cascade fea- 
tures in the spectrum, we set the required signal-to-noise ratio 
to be 5cr per logarithmic energy bin width of AE/E = 0.1. 
This corresponds to a 19er detection for integrated flux. From 
the expected flux in our standard model for M82, the required 

1 CTA: http://www.cta-observatory.org/ 
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FIG. 3. — Left panel: Gamma-ray spectrum of M82. The solid curve is the total flux, i.e. absorbed plus cascade component. Three model predictions for the 
intrinsic (no absorption), absorbed, and cascade components are also shown by the curve marking shown in the figure. We adopt (F, T„ max ) = (2.4,512 TeV) 
for the intrinsic spectrum. For i nterstellar radiation field SED model, we adopt the SK07 model. The observed data of Fermi I Abdo et al. 2010a) and VERITAS 
( VERITAS Collaboration et al. 2009) are also shown. Right panel: Same as the left panel, but for NGC 253 with (T,T Pim!lx ) = (2.3,512 TeV). The observed 
data of Fermi I Abdo et al. 2010a) is also shown. The single VHE flux point for NGC 253 is computed from the integral photon flux >220 GeV reported by the 
H.E.S.S. ( Acero et al. 2009) assuming a power-law spectral model with photon index between 2.0 and 3.0. 
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FIG. 4. — Same as Fig. [3] but with the total flux for various T, 7p, max , and interstellar radiation SED models. Left panel: Solid, dashed, dotted, and dot-dashed 
curve is the model for (T,T Pimat ) = (2.4,512 TeV), (2.3,512 TeV), (2.4,45.3 TeV), and (2.4,4.0 TeV) with the SK07 interstellar SED model, respectively. Double 
dot-dashed curve is the model for (r, 7},, ra ax) = (2.4, 512 TeV) with the dCdP09 interstellar SED model. Right panel: Solid, dashed, dotted, and dot-dashed curve 
is the model for (r,r p , max ) = (2.3,512 TeV), (2.2,512 TeV), (2.3,45.3 TeV), and (2.2,4.0 TeV) with the SK07 SED model, respectively. Double dot-dashed 
curve is the model for (T,T PiaaL ) = (2.3,512 TeV) with the DST05 SED model. 



observing time to achieve this signal-to-noise ratio is 14, 19, 
and 330 hours for 300 GeV, 1 TeV, and 3 TeV, respectively. 
In the case of NGC 253, the required time is 14, 14, and 180 
hours for the above 3 energy bands in our standard model, re- 
spectively. Thus, it would be possible for CTA to obtain high 
energy resolution spectra to see the cascade effects within rea- 
sonable observational times. 

5. CONCLUSIONS 

In this paper, we estimated the contribution from the cas- 
cade emissions which are resulted from the inverse Compton 
scattering of the interstellar radiation field in the galaxies by 
the e + e~ pairs created by 7-ray attenuation. To evaluate the 
7-ray attenuation, we first modeled the 7-7 optical depth of 
very high energy 7-ray (>30 GeV) absorption by the interstel- 
lar radiation field via positron and electron, e + e~, pair creation 
in the two nearby bright starburst galaxies, M82 and NGC 
253. To take account of the interstellar radiation, we adopted 
th e SED model in the cent r al star forming region provided 
by Siebenmorgen & Kriigel (2007). Since 7-rays are created 



inside of the central star formation region, absorption effects 
from such obscured emission is important. We found that the 
attenuation effects becomes significant above ~ 5 TeV and 
~ 9 TeV for M82 and NGC 253, respectively. 

By using the interste l lar ra diation field SED model by 
Siebenmorgen & Kriigel (2007) and our 7-7 optical depth 
model based on it, we found that the cascade emissions would 
be a probe of the existence of high energy cosmic -rays (> 10 
TeV) in starburst galaxies. The resulting total flux including 
cascade emission would be ^18 % and ~45% higher than the 
absorbed flux above 1 TeV with the maximum proton kinetic 
energy, 45.3 TeV and 512 TeV, resp ectively. Even when we 
use the Ide Cead el Pozo et al. (2009) in terstellar SEP model 
for M82 and the Domingo-Santamana & Torres (2005) inter- 
stellar SED model for NGC 253, the effect of cascade emis- 
sion is not significantly changed. These differences would be 
a new indirect investigator of the highest cosmic -ray energy, 
although it would be difficult to see the direct 7-ray signature 
because of the attenuation. It might be possible to investigate 
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this cascade signature through a detailed differential spectrum 
observation with future high signal-to-noise ratio 7-ray obser- 
vations such as CTA by comparing with the theoretical 7-ray 
emission models including this cascade effect. 
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